Weak asymptotic of shock wave formation process 
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Abstract 

We construct an asymptotic (in a weak sense) solution corresponding 
to the shock wave formation in a special situation. 

1 Introduction 

We consider the problem of shock wave formation for the following Hopf type 
equation: 

where we assume that / e C'^ and the inequality f"{u) > holds on the range 
of the solution u. We shall consider the special initial condition for Eq. (1): 

^t=o ^ "0 + - u°a)H{ai ~ x) + {U - ui{x))H{a2 - x), (2) 

where Wq, and ai > 02 are constants, H is the Heaviside function, the 
function ui{x) is determined by the equation 

f{ui{x)) = -Kx + h, K,b = const, (3) 

and, in addition, we assume that wi(ai) = Uq and Ui(a2) = U. 

Such a function appears^ in the construction of the entropy solution to the 
Cauchy problem with an "unstable" initial jump. 

It follows from the choice of such an initial condition that the approximation 
of problem (l)-(2) (a weak asymptotic solution) for all t is an element of the 
asymptotic subalgebra 

B{1, Hi{x - V5i, e), H2{x - v?2, e)} 

introduced in 0. 

Roughly speaking, this means that at any time moment the weak limit of 
the weak asymptotic solution is a linear combination of the Heaviside functions 
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H{x — (fii) and H{x — (^2) with smooth in t (e > 0) cocfBcients and that there 
are no additional jumps. In turn, this means that at time 



t* 



ai - a2 
f'iU) - /'(wO 



all characteristics meet at the same point x* = ai + Vit* , Vi = ('{uq), V2 = 
f'{U),i = 1,2. 

More precisely, for < f < t*, the solution of problem (l)-(2) is given by 
the formula 



u = Uq + (ui{xo{x,t)) - u*^)H{tpi - x) + (U - ui{xo{x,t)))H{tp2 - x), 
where the function ui(xo(x,t)) has the form 



(4) 



ux{xQ{x,ty) = Mil ai + 



V'o 



Ml 



X - bt 



Here (pi{t) = + Vit, i = 1,2, -00 = Vi{t) — 'P2{t), and t/jq — ai — a2. 
For t = t* the plot of the function u = u{x, t) is the graph 

((-00, x*), U) U {x*, {U, u°)) U ((a;*, 00), u°). 

We note that if we set 



' ui{xo{x,t)), t<t*, t>t*, 



u{x,t) 



^01 



t = t*, x<x*, 
U, t = t*, x>x*, 

Ji e K, U], t = t*, x = X*, 



then the function u{x, t) is defined for all values of t and, for t <t*, is a solution 
of Eq. (1) satisfying the initial condition (2) for t = 0. Our goal is to "correct" 

the function u(x, t) and to obtain an analytic formula that, for t < f* , determines 
a function close to u{x,t) and, for t > t* , a function close to the function 



u = ul + H{c{t - t*) - (x - x*))U, 



(5) 



where 



[/(")] 



x=ct 



f{u)-f{<) 
u-4 ■ 



The answer is given by formula (10) below. 

We note that the function u determined by relation (4) for t < f* is contin- 
uous everywhere except the points lying on the curves x = ^Pi(t), i ~ 1,2, and, 
at points of these curves, the function has weak discontinuities (the derivatives 
of the function have jumps at these points). Therefore, the formation of the 
shock wave (5) from (4) can be treated as the result of interaction (confluence) 
of weak discontinuities. Moreover, for t < t* , although the derivatives are dis- 
continuous, the solution of such problems (that is continuous, but with jumps 
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of the derivatives on some smooth nonintersecting curves) can be constructed 
by the method of characteristics. 

We note that, in the case f{u) = , the problem of constructing the global 
asymptotic solution of problem (1), (2) was solved as an example in 2|. The 
asymptotic solution constructed in (2 is a weak asymptotic solution. We recall 
how it is determined. By Ot)'{s") we denote generalized functions that, in 
general, depend on the parameters t and e and are such that for any test function 
rj{x), the estimate 

(Oi,,(£"),r;(x))=0(e") 

holds, where the estimate on the right-hand side is understood in the usual sense 
and locally uniform in t, i.e., |0(e°)| < Cte" for t e [0,T]. 

A function = u^{x, t) is called a weak asymptotic solution of problem (1), 
(2) if 



dt dx 



= Ov'{e). (6) 



The goal of this paper is to construct such a function in the case of a general 
convex nonlinearity f{u). This is achieved in Sections 2 and 3. 

In Section 4 we introduce auxiliary formulas and statements of the weak 
asymptotic method. 

We note that if the solution satisfies the Oleinik-Kruzhkov stability con- 
ditions [3 ^ , then it follows from (6) that differs from u by a measure |S] 
whose values are estimated as 0{e). Indeed, it is easy to verify that the right- 
hand sides in (6) arising in our construction belong to C([0, T], L^(iRj)) and 
can be estimated as 0{e) in the sense of the L^-norm. Therefore, according to 
the results in |31 0] , ite is an asymptotic of the solution to the Cauchy problem 
(l)-(2) in L^. This is done in Section 5. 

We also note that the asymptotic (in the usual weak sense) solution de- 
scribing the global behavior of the solution of the Cauchy problem with a small 
viscosity and a smooth initial condition for the equation 

du df{u) d^u 
dt dx dx^ 

was first constructed by A. M. Il'in 6 . This was an important achievement in 
the asymptotic theory. 

In contrast to our paper, in A. M. Il'in's paper an arbitrary smooth initial 
condition was considered. In in the case f{u) = m^, it was explained how 
the solution constructed there can be used to obtain the global weak asymptotic 
for a more general Cauchy problem. For this, it was proposed to consider an 
interpolation of the initial function by linear splines. Also, we can use different 
approach. For given smooth initial data uq{x), a; G R, we can find (assume 
finite) set of points a;° G R, A; € 1, ---jN, N £ N, which reaches the point of 
the gradient catastrophe in the moment t — tk- Now, instead of given initial 
data uo{x), x G R, we impose initial data uoeix), x G R, which differs from the 
function uo{x), a; G R, in the intervals [ e''], 0< ^i<l,ke 1,...,N. 
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In those intervals the function uos{x), a; G R, has form (10). It is obvious that 
we have 

||uo(a;) - uoe(a;)||Li(R) = 0{s^). 

Then, we solve the Cauchy problem corresponding to new initial data uoe{x), 
X £ R, using method of characteristics in a way that for the " inserted" parts (the 
one in the intervals [x^ — e'^ ,x^ + e'^], < /j. < 1, k ^ 1, N) we use "the new 
characteristics" given by (11) and for the rest of the function uoe{x) ( = uo{x) 
for X ^ [x^ — e^.Xfe + ^^]) use ordinary characteristics. This will be the 
subject of further investigations. 



2 Description of the formula 

for the weak asymptotic solution 

To construct a weak asymptotic solution describing the passage from (4) to (5), 
we introduce some auxiliary constructions. 

We define a function ^{xq) as a solution of the implicit equation 

U + u°o = ui{xq)+ui{^{xo)), (7) 

which is solvable due to (3). 

Obviously, ^ : [02, fli] [02, fli] is a smooth isomorphism and ^(^(xo)) = Xq. 
We introduce the function Ui (xq , p) , by setting 

Ui{xo, p) = B2{p)ui{xo) + Bi{p)ui{^{xo)), (8) 

where the functions Bi{p), i = 1,2, are defined in Lemma 4.1 and the function 
p = p{t) is defined below, see (12), (13) and 

T= , ipio{t) = ai + f {ui{ai))t, 1 = 1,2. 

Note that, by (7), (8), and the formulas for Hi at the end of Lemma 4.1, we 
have 

U + ul-Ui{xQ,p) = Ui{^{xq),p), Ui{xo,p)=ui{xo)+O{p~^), p ^ 00. 

(9) 

We shall seek a weak asymptotic solution of problem (1) (2) in the form 

Ue{x,t) =Uo+ {y\{xo{x,t,T),p) - Uq^UJi( ^^ ^ ^ 

+ {u-U^{x^{x,t,T),p)y2{^^y (10) 

where uJi{z) — > 0, 1 as 2; — > =F00j ^ ^ = 0{\t\~^), where j^;] — > 00, a > 

and N > are arbitrary numbers, and = (j)i{t,e), i = 1,2, xo{x,t,T) are the 
desired functions. 
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As noted in the last Section (see Sec. 4.1), the functions uji{{4'i — x)/e) 
approximate (in the weak sense) the Heaviside functions H((j)i — x), 

= H{(t)i-x) + Ov'{£), i = l,2. 



£ 

We shall seek the functions (t)i — (piit, e), i = 1,2, in the form 

(j)i=ipi{t,T)+ll}Q^{t,T), i=l,2,. 

Here, (f){t, t) is such that it satisfies (j){t, 't)\_^^^ = 0. Furthermore, '^i{t, t) is an 
analog of the trajectories of weak discontinuities of ifiit) in (4). The functions 
(Pi{t,T), i = 1,2, can be found from the equations for the "new characteristics" 
and, as T ^ oo (i.e., before the confluence of weak singularities), these functions 
are close to the trajectories (pi{t) from the preceding section. 

To find the functions Xo{x,t,T), we introduce the differential equation for 
the "new characteristics" 



— = B2{p)fmxo,p))+B,{p)f'{U^{axo),p))+q{T,p), 



= xo. (11) 



t=0 

The function q{T, p) is assumed to be smooth and to satisfy the estimate 

|Tg(T,p)| < const. (12) 

Its appearance itself is caused by the fact that the function Ui{x,p), which 
replaces the function ui{xq) in formula (4), depends on time (via the function 
p = p{t) determined in (13), (14)). Therefore, this function is not preserved 
along the usual trajectories corresponding to quasilinear equations. The "new 
trajectories" are just given by Eq. (11), where the function p is determined as 
follows. By x{xo,t, r) we denote the solution of (11) and introduce the functions 
(pi{t,T) = x{ai,t,T), i = 1,2. We set 

^ 0l{t,T) - ip2{t,T) ^ <^l(^,r) - (f)2{t,T) 
£ £ 

We note that Ui{ai,T) = B2U^ + BiU and Ui{a2,T) = B2U + BiUq; hence 
from (11) we easily obtain the following equation for p = p{t): 

g = {B2{p) - B,{p)){f'{B24 + BiU) - f'{B2U + B,u',)){i,',)-\ (13) 
Obviously, by definition, 

pT~^ 1 as T — » 00 and 3^ > (since Vo < 0)- (14) 

(IT 

We denote the right-hand side of (13) by G{p). Obviously, G{po) = 0, where po 
is a number such that i?i(po) = -82(^0)) and hence (see Lemma 4.1) 

Bi{po) = B2{po) = 1/2. (15) 
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We assume that po > 0. This is a condition imposed on Bj. It is easy to verify 
= 0, while 



, dG 
that — 

dp 



P=PO 



fG 
dp-^ 



P=PO 



= -8B^^,(f/-«°)/"(^^)^0. (16) 



It follows from Eq. (13) and inequality (16) that the relations 
p^po + 0{l/\r\), p = 0{l/\Tf) 



(17) 



hold as T ^ —oo. 

Thus, independently of (11), the function p = p(t) is defined as a solution of 
problem (13), (14). Therefore, the function x{xQ,t,T) from (11) is also defined. 
Now put 

x{xo,t,T) = X{xo,t) + V'o-'^i(a;o,T), (18) 



where 

X(xo, f) = xo + f'{uoixo))t = a;of/'o(V'o)~^ 
Inserting x in (11) instead of x we have 



bt. 



=— r 

V'o'T Jo 



B2{p)f'{U^{xo,p))+B,{p)f'{U^{^{xo),p)) 



+ q{T',p)-f'{uo{xo)) 



dr' 



— f 

%T- Jo 



B2{p)nUi{xo,p))+B,{p)f'{U,{^{xo),p)) 
+ g(T',p)]rfr'-(Vg)-^xo-6(V^)-^ 
It is easy to verify that the following representation is true: 



(19) 



X = X + 



^ r 

''or Jo 



B2{p)nUi{xo,p))+B^{p)f'{U^iaxo),p)) + q{T',p)\ dr', 
which follows from the identity 

xoMi'o)-' + bt- M^o)~'xo - V-oKV-o)"' 



-6^ = x*. 



It is not difficult to see that the solution x given by formula (19) is not the exact 
solution of (11). Actually, for t = (i.e. for r — > +oo)from (19) we obtain 

{Xo + yJoXi)\t=o^xo + 0{e). 

Obviously, for f G [0, T], T G R, we have: 

x{xo, t, t) = x{xo,t, t) + 0{s). 
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It is easy to verify that the term 0{e) in the last relation has the form 



0(£) = Vo^l 



t=0 ijj'oT Jo 



B2 (p) /' (C/i (xo , p) ) + (p) /' (C/i (e(xo ) , p) ) 



+ q{t', p) - f'{uo{xo)) dr' =^ sg{xo) 



Finally, we obtain 



x(xo, t, t) = x{xo,t, t) + egixo). 
Let us calculate the derivative By (18), (19), we have 



dx _ -00 r d 
dxo %T Jo dxo 



B2ip)nUi{xo,p)) + B,ip)fiUiiaxo),p)) 



dr 



i'oT- Jo 



,1 dUi 



B2{p)f"{U^{xo,p)) - B,{p)f"{U^{^{xo),p))\ -^{xo,p)dT 

(20) 



Here we used the relation 

^^{xo,p) = -^^i^ixo),p), 

OXo OXq 

which follows from the definition of the function Ui{xo,p) in (9). 

We agree that the symbol ~ denotes the following equivalence relation 



/ 

f 9 ^ lim — = const ^ 0. 



(21) 



Then, as r ^ — oo, we have 
dUi 



dxo ^ 2 r' 



Bo- 



1 1 



2 r' 



U^{X0,p)^^^;^, U,{aXo),p) 



Therefore, 

'B2{p)f"{Ui{xo,p)) - Bi{p)f"{Ui{axo),p)) 
Hence the integral in (20) converges as r — > — oo and 



1 

T 



dx 

dxo Ip'oT ' 



T — »■ — oo. 



As T — »■ 00, we have Ui{xo,p) ui{xo) (since B2 — » 1) and the integrand 
in (20) tends to the limit 



Wo 



dui 
dxo 
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00. 



Thus 

dx Vo 
dxo ip^ ' 

Wc note that the solvabihty of the equation x{xq, t,x) = x with respect to xq 
globally in t can hardly be ignored. 

In our constructions, we shall hence use the following approximate expression 
for the solution of Eq. (11), namely, 

x{xQ,t, t) = x{xo, t, t) + £{g{xo) + Axq), 

where A> Q, A = const. Clearly, we have 

— Xq+ eAxq, 

t=0 

and hence xo{x,t,T)\t=o = x — sAx + 0{e). 

As is easy to verify, this means that, in the sense of Oi>'-estimates, the initial 
condition in (???) will be satisfied with accuracy up to Ot>'{e)- We prove that 
the constant A can be chosen so that the inequality 



dx 
dxo 



> 



holds uniformly in t. 
We have 



OXo %T Jo 



B^f"{Ur{xo,p)) - Brf"{Ur{^{xo),p)) 



Ou 1 

-/"(ui(a;o))^J dr' +eg'{xo)+sA. 

Recall that t* = k~^, ipo{t*) = 0, r = il)o{t)/e. Hence for t < t*, by Lemma 4.2, 
we have the estimate 



1-kt- 



B^f"{Uiixo,p)) - B,f"{U^{^{xo),p)) - f"{u,{xo))p-] dr' = 0{e) 



^'oT JO 

Similarly, for t >t*, we have 

^0 



%T Jo 



B2f"{Ui{xo,p))-B,f"{Ui{^{xo),p))-f"{ui{xo))g^^ dr' = 0{s). 



It follows from these estimates that there is a possibility to choose the con- 
stant A. Thus the equation 



Xo{xo,t) + ipoXi{xo, t,s) + e{g{xo) + Axo) = x 



(22) 



can be globally solved with respect to xq. 

In this case, the derivatives of the exact solution of Eq. (11) differ from the 
function in the right-hand side of (22) and from the function x{xo, t, r) by 0(e). 
Therefore, in what follows, to simplify the calculations, we shall use all these 
functions. 
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3 Construction 

of the weak asymptotic solution 

Wc substitute the function Ug{x,t) into Eq. (1). Using Lemma 4.1 and the 
formula for weak asymptotic of the approximations in Sec. 4, we obtain 



du, d 

-dt + a^^^"^ 



(j)lt{Ul{XQ{(t)l,t,T),p) - Uq)5{x - ^i) 

+ 4>2t{U - Ui{xo{(f)2,t,T),p))6(x - (j)2) 

dUi dxo 



dxo dt 
dU\ dxo 



H{4>^ -x)- H{ct>2 - x)] 
[B2{p)f'{Ui{xo{x,t,r),p)) 



+ 

dxQ dx 

+ B,{p)f'{U,{^{xo{x,t,T)),p))] [H{^, -x)- H{<j>2 - x)] 

- 5{x - ct>r) [b2{p) (/([/i(a;o(<Ai, i,T),p)) - 

+ i3i (p) (/([/)- /(C/i (e(xo (01 , t, r)), p))) " 

- 5{x - 02) [Bi{p) [f{UMMh:t, r)), p)) - f{u°,)) 
B2{p){fiU) - f{U^{xoih,t,r),p)))' 

[H{(l>^-x)-H{ct>2-x)]+Ov'{e). 

(23) 



+ 



+ ^(a=o,p) 



Xo=Xo{x,t,T) 



Although there are rather many terms on the right-hand side, it is easy 
to understand this formula. The terms containing the factors {H{<pi — x) — 
H{4>2 — x)) correspond to the substitution into the equation of the function 
determined in (10) between the points x = ipi and with Lemma 4.1 taken into 
account. 

The terms containing the delta functions, i.e., the factors 5{x — 0i), i = 1,2, 
appear due to the fact that U{xo{(j)i{t,T), p)) ^ Uq and U{xo{(j)i{t,T), p)) ^ U, 
but as T — > oo (i.e., before the interaction) we have p ~ r (see (14)) and hence, 
for any A'' > 0, we have 

U{xa{Mt,r),p))-ul = 0(£^), U{xo{cl>2{t,T), p)) -U = O(e^). 

We start analyzing the terms in (22) from the last one (which has the estimate 
0{e~^) in the C-norm): 



= e-^i;'op[B'^pU,{xo) + B[pU,{C{xo))] 



= e-^i^'^pB'^p [U + ul~ 2ui {xa{x, t, r))] . (24) 
Applying the Taylor formula at the points x = 4>\ and x = (j)2, for any test 
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function r]{x), we obtain 



J d)2 



= YPpS2p(C^ + ^^o) -<f>i)+5{x- <f>2)}, v{x)) + 0(e), (25) 

and = 0{\p\~^) for any > as p ^ oo, i?2p const as p ^ po 
(r — > —do), and p = 0(|r|~^) as r ^ — oo (see (17)). 
Let us consider the remaining term. We have 



J d>2 



- dx 

ui{xoix,t,T))ri{x) dx = / ui{xo)r]{x{xo,t,T) + ^o(p)- — dxo. (26) 

Ja2 "'^0 



Let us note that (see (18)) 



^X f/iO\-l 



Hence the right-hand side in (23) is bounded in the weak sense as e — > 0. We 
now note that the following relations hold: 

r]{x{xo,t,T) +ipo<l>) = v{xiO'i,t,T) + 'tpo<i>) + v'xT— 

OXo 

Ci e (ai,a;o), i = 1,2. 

Recalling that x{ai, t, r) = (pi{t, r), i = 1,2, and again using (26), (14) and (17) 
we obtain 

s~^'tpQpB'2 / ui{xo{x,t,T))r]{x) dx 

J <t>2 

= ^ {d{x - + 5{x - (l)2),v{x)) 



X £ ui{xo) (^(^0°)"' + dxo + 0{e). 



Finally, we have 

dUi 



at 

where 



[i?(0i -x)- H{4>2 - x)] = g{T,p){6{x -<^i) + 6{x- <^2)) + 0-D'{e), 



9{r, P) = ^PPB'2,{U + u'o) - <PrB'2, £ m{xo) ((^g)"' + ^^^o- (28) 

It is easy to see that, by formula (17), we have the estimate 

\t'^9{t,p)\ < const. 
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Moreover, the function g{T,p) is integrablc, and the integral g{T, p) dr con- 
verges. Indeed, the integral of the first term converges because of the estimates 
given after formula (24), and the integral of the second term, in its properties, 
coincides with the last integral in formula (20). 

Now we consider the remaining terms that contain the difference H[(j)i — 
x) — H{(f)2 — x) as the multiplier. For any function ri{x) G Cq^, taking into 
account the relation 

dxo dxQ dx 
dt dx dt ' 

we have 



dU\ dxo dUi dxo 



dxo dt dxo dx 

+ B[{p)f'(U + u^o-Ui{xo{x,t,T),p)) 



(B2{p)f'{U^{xo{x,t,T),p)) 

{H{(t>i-x)-H{ct>2-x)),ii{x) 



dUi dxo 



02 dxo dx 
+ Bi(p)/'(C/i(C(a;o(x,i,r)),p)) 



-^+B2{p)fmxo{x,t,T),p)) 



T}{x) dx. 



(29) 



By (11), the expression in square brackets on the right-hand side of (28) is 
just q{T,p). 

We consider the integral 

dUi dxo 



I 

J (h-2 



-r){x) dx 



dxQ dx 

and pass to the variables xq precisely as in (25). We obtain 



/ 

J(h 



dUi dxp 
, dxo dx 



r]{x) dx 



dUi 



r]{x{xo,t, t) + tpoi') dxo- 



Recall that 



dUi 
dxo 



1 

T 



-DO, 



a2 

dU\ du\{xQ) 



00. 



dxo dxo 

Prom the conjectural estimate (12) for the function q{T,p), using the Taylor 
formula as in (24), we obtain 

dUi dxo 



02 dxo dx ' 2 



dxo {5{x - <l)i) + 5{x - (j)2), ri{x)) + 0{e). 

Taking into account the definition of the function U\{xo,p), we can easily 
calculate the integral on the right-hand side of the last formula and obtain 



-/ 

J(h 



dUi dxo 



q{B2 - Bi)(wO - U) 



{S{x - ct>i) + 5{x - ci>2), v{x)) + 0{e). (30) 
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We choose the function g(r, p) so that the following relation hold: 



q{B2-B^){ul-U) = -g{T,p). (31) 

Obviously, we have 

g(T,p)^.9(r,p) = 0(r-^) ViV, r ^ oo; 

q{T, p) Tg{T, p) r-. - , r ^ -00. 
r 

Hence the estimate (12) holds and our constructions that lead to (29) are well 
defined. 

It is left to obtain the function <f) appearing in the definition of the functions 
phii, i = 1,2. To do that we will use the results from Section 4.3. Equating 
with zero the remaining coefficients oi S{x ~ (pi), i ~ 1,2, (only such expressions 
(mod 0'pi{e)) are left on the right-hand side of (22)), we obtain 

<Pit (;7(a:o ((/.I , i , r) , p) - yO) - B2 (p) (/ (C/i (.To {<f>i,t,T), p)) - f{u°)) 

-i?i(p)(/(C/)-/(C/i(e(xo(</.i,t,r)),p))) =0, (32) 

ht{U - C/i(xo(02, t, t),p)) - B2{p) {f(U) - /(C/i(xo(02, t, t), p))) 

-Bi(p)(/(C/i(^(To(02,t,T)),p))-/(w[!)) =0. (33) 

According to (47) we have to prove that preceding equations are correct when 
T — > ±00 and to find (p such that their sum be equal to zero. 

By the definition of the functions <pi{t,T), i = 1,2, as t oc (i.e., before 
the interaction), the limit of the expressions on the left-hand side of relations 
(29), (30) is equal to zero, and these relations admit the estimate 0{t~^) for 
any A/' > as r — > 00. This follows from the relations: p/r — > 1 as r — > 00 and 

B2 = l + 0(p-^), Bi=0(p-^), p^oo (r^oo). (34) 

We write the limit of these relations for r — > —00. Recall that 

Bi{p) = \+0{\T\-^), p = po + 0(|r|-i), r^-00, i = l,2. (35) 
Therefore, denoting the limit of (pit as r — > —00 by i = 1, 2, we obtain 
<t>u{^^^)=\i.f{U)-f{ul)), i = l,2, (36) 

or 

Denoting, as usual, 

f{U) f{ul) ^ [/] 

u-ul M' 
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wc can determine the general limit (j) {t) of the functions (j)i{T,t), i = 1,2, as 
r ^ — (X) by the relation 

r = r(t*) + |{|(i-n- (38) 

Relations (36) (or (38)) mean that, for t > t* , the trajectories x = (pi and 
X = (p2 are close to the line 

i.e., to the trajectory of the shock wave (5). 

Let us investigate the trajectories x = (j>i, i = 1,2, in more detail. 

By Lo{z) we denote the function satisfying the same conditions as the func- 
tions coi, i = 1,2, in (10). 

We prove that the following relations hold: 

(f>i{t,T)-4>i{t,T) = 0{s), i = l,2, (39) 

where 

Mt,r) = (l-co(r))93,(r,t)+c.(r) ( x* + ^/]{t - t*) 

Here <^i(r,i) = X{ai,t) + ■ipoXi{ai,T), (see (18)), x* = ipio{t*) = ^p2o{t*), 
and (f)i{t,T), i = 1,2, are the desired trajectories of singularities determined 
(mod 0{e)) by the relations 

(j)i{t,T) = 0i{T,t) +tpo4>- 

To prove (39), it suffices to set (p^it*) = x* in (38) and to note that 
(pi{0,t*) = X* , i = 1,2. It remains to note that the functions (pi{T,t) can 
be represented in the form 

^i(r, t) = x*+ Mi^'or)-' [[B^nUi^a,, p)) + Bif{U{i{a,),p)) + p)] dr' . 

Jo 

(40) 

This follows from (18), (19) with the relation x* = bt* = — 6Vo/V'o taken into 
account (see formula (3)). 

Now we apply Lemma 4.2 and see that relation (39) is proved. The state- 
ment wc have proved means that <pi from (39) provide a family of expressions 
for trajectories close to those trajectories we want to construct. These approx- 
imate trajectories, with accuracy 0{s), are independent of the choice of the 
function uj(t). It is only required that this function satisfy same conditions as 
the functions u>i, i = 1,2, from (10). 

Now let us calculate the function xo{(f>i,t,T), i = 1,2. By definition, this 
is the initial point of the trajectory x = (j)i{t,T), i = 1,2. Clearly, for t < t* , 
we have (j)i{t,T) = (pi{t,T) + 0{e) and xo{iPi,t,T) = ai. For t > t* , we have 
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(j)t{t, t) - (j) (i) — > as e ^ 0. By relation (37), for ^ (t*) = x*, we see that in 
this case the initial point is 

x = (p-{0) = x* -^/]t*. 

[u\ 

By the inequalities f'{U) < [f]/[u] < /'(uq), this imphes that x G (02, ai). 
We set 

XQ{4>i,T) =ai + n{T){x-ai), «=1,2, 

where ri(T) is some (generating) function satisfying same conditions as the func- 
tions cui, i = 1,2, from (10). 
Let us prove the relations 

Mt, r) - {x{Xo{cj)i,T),t, t) + ^o4>) = 0{e), i = 1,2. (41) 

We restrict ourselves only to the case i = 1. We have 

Ui{ai,p) = Ui{Xo,p) - 0(i: - "'i)^"!^! + "^(* ~ ai),p). 



Ui{^{ai),p) = U + u°-Ui{aup) 

) — ('T. — rr..- ^, 

dxo 



= U + Uq- Ui{Xo,p) -{x- ai)- — {ai+aD.{x - ai),p). 



where a e (0, 1). 

From these relations, formula (40), and representation for x from Section 2, 
we obtain 

0^{t,T) - x{XQ{<l>x,T),t,T) (42) 



Jo 



S2^(C/i(ai,p);C/i(Xo,p)) 

Bi^{U + u°- Ui{ai,p); U + u° - Ui{Xo,p)) 
du 



where 



dU 

X f2— — (ui + aQ.{x — ai),p) dr' 
0x0 



5u ' ' A-B A^B ■' ^ ' 



We note that the integral on the right-hand side of (42) converges as r ^ +00 
because the function Q is contained in the integrand. The convergence of the 
integral as r — > — 00 can be verified in the same way as the convergence of the 
last integral on the right-hand side of (20). Hence, by Lemma 4.2, we have 

if lit, t) - x{Xo{4>i,T),t,T) = 0(e), 

and hence, by (38) and (39), we obtain (41). From (41) we obtain the relation 

Ux{xo{(l)i,t,T),p)-Ux{Xo{<j)i,T),p)=0{s), i = l,2. (43) 
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By construction, the limits of the expressions on the left-hand sides in (31) 
and (32) are equal to zero as r ^ cxd (i.e., before the interaction). Moreover, the 
difference between the limit and the prelimit expression is 0{p~^) = 0(t~^) 
for any > 0. 

By (32), these expressions also tend to zero as r ^ —oo, and the difference 
between the limit and the prelimit expression is 0{Bi — 1/2) = 0{p — po) = 
0{\t\~^), t — oo. Therefore, by the results of Sec. 4.2 about the linear inde- 
pendence, for the sum of terms with (5-functions in (22) to admit the estimate 
0-D'{e), it is sufficient that the sum of expressions on the left-hand sides of (31) 
and (32) be equal to zero. Thus we obtain the equation 

= B2{p){f{U,^,)) - f{ul)) + B,{p){m - /(C/id))) 

+ B2{p){f{U) - /(t/i(2))) +Bi(p)(/(?7i(2)) - M)). (44) 

Here, for brevity, we denote 

Ui{i) = Ui{xo{(pt,t,T),p), U^i) = Ui{£,{xo{(j)i,t,T)),p), i = l,2. 



We note that 



ifit +ipo^{T(p), i=l,2. 



We agree to denote f k, g \i 



lim - = 1. 
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It is easy to verify that as t oo, we have 

U - C/i(2) « t^i(i) -ul^U- Z7i(i) « Z7i(2) - « Bi(C/ - «°). (45) 
Similarly, 

/(f/id)) - f{ut) « nul)B,{U - ul), (46) 

/([/)- /(?7i(2))«/'(f/)Bi(f/-«g), 

/([/)- /(C/i(i))«/'(C/)i3i(t/-n°), 
f{Uri2))-f{ul)^f{ul)B^{U-ul). 

Next, by (16), we have B'2{p) ~ 1 — B2 and hence the relation 5 ~ 1 — B2 
holds as T ^ +00. 

As r — > —00, the coefficient of -^{Tcj)) in Eq. (44) is equal to U — Uq ^ 0. 
Therefore, Eq. (44) is solvable for ^ and its solution is a bounded function 

decreasing as t ^ 00. 

To write the solution of Eq. (44), we note that, with accuracy 0(e), by (41), 
we can replace the arguments xo((j)i,t,T) by XQ{(f)i,T) in the functions CAq), 
and by (38), the function Xo(^j,T) can be determined actually independent of 
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the functions (pi (everywhere here i,j = 1,2). Hence Eq. (44) is indeed a linear 
equation with respect to cj) and its solution can be easily found. 
This solution has the form 

^ = ii^'or)-' r {U-u°-U,{Xo ih ,t),p) + U, {Xq (</>!, t) , p)) 
Jo 

X ( - <^2 [f/ - C/i (Xo (</.2 , r) , p)] - ^1 [f/i (Xo (</.! , r) , p) - 

+ {B2(p)(/(C/i(i)) - M)) + B^{p){f{U) - /(C/id))) 

+ B2{p){f{U) - /([/i(2))) + Bi(p)(/(?7i(2)) - f{ul))}) dr'. 

By (45) and (46), the integral on the right-hand side in the last relation converges 
as T ^ cx) and ^ = 0{t~^) as r ^ oo. 

4 Auxiliary formulas and statements of weak 
asymptotic method 

4.1 Nonlinear superposition of approximations of Heavi- 
side functions 



Suppose that (Jj{z) — > 0, 1 as 2; ^ —00, 2 — > 00, ^ J = 0{\z\ ), j = 1,2, 

\z\ oc, N is a sufiicicintly large number, and tpi, 1^2 are some continuous 
functions of the variable t. 

It is easy to verify that the functions u)j{{x — (j){t))/e) approximate in the 
weak sense the Hcavisidc function H{x — (pit)). Indeed, the properties of the 
functions LOj{z) imply the relations 

u,j{z)-H{z) = 0{\z\-^), N>0, j = l,2. 

Hence, for any test function ip{x), we have 

' .T — (6 ^ 

£ 

Lemma 4.1. For any - function f(x), the following relation holds: 



H{x-ip),i)^ =s J {u}j{z) - H{z))'ip{(p + sz) dz = 0{e). 



/(a + 6a;,(^j+6a;2(^ 

= /(a) + H{ip^ - x){B2{f{a + b)- /(a)) + Bi(/(a + b + c)-f{a + c))} 
+ H{ip2 - x){B,{f{a + c) - /(a)) + B2{f{a + b + c)-f{a + b))} + Ot,,{s), 

where 

Bj=Bj( j^~'^^ y J = 1,2 Bi+B2 = l, 
B2{z) 1 as z — > cxD, 52(2) as z ^ —00. 
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Proof: First, we prove the relation 



/ a+buJi 



'CUJ2 



tp2- X 



= f{a+bH{ipi~x)+cH{^2~x))+Ov'ie). 



Indeed, we have 

/(a + buji + CUJ2) = f(a + bH{ipi - x) + c{ip2 - x)^ 

+ /' (a + - H{^i - x))) + c{u;2 - H{^2 - x))^ 

X [{uji - H{ipi - x))b + {uj2 - H{ip2 - x))c]. 

Now we verify that if 17(2;, (p, e) is a bounded function, then 

g{x,ip,e)[uji - H{(pi - x)] = Ov'{e). 

For any test function ipix), we have 

f (pi ~ X 



9ix,(p,s) 



ip{x) dx 



H{ipi - x) 

£ J 9i'^i+£z,ip,e)[uj{z)-H{z)]ij{ip-ez)dz 
< econst / \lu{z) — H{z)\ dz. 



This imphes 

/(a + buji + /(a + bH{ipi - x) + cH{ip2 - x)) + Ov (e). 

Next, it is easy to verify the relation 

f{a + bH^ + CH2) = f{a) + H,[f{a + b) - f{a)] + H2[f{a + c) - /(a)] 

+ H1H2 if {a + + c) - /(a + c) - /(a + 6) + /(a)) , 

h/^' H{ipj-x), J -1,2. 
It remains to note that we have 

H{ipi - x)H{if>2 -x)^ BiH{ipi -x) + B2H{ip2 - x) + Ov'{e), 



(l>l{z)uJ2 Z 



dz. 



B2 ^l-Bi. 



For the proof of these and similar relations, see m21[7]. The proof of the lemma 
is complete. □ 
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4.2 Asymptotic linear independence 

If we want to consider linear combinations of generalized functions with accuracy 
O-D'i^"), then we need to modify the notion of linear independence. This mod- 
ification plays the key role in considerations related to the soliton interaction 
problem. 

Indeed, let <j)i ^ 4>2 be independent of x. We consider the relation 
gi5{x - 0i) +52<5(x - (j}2) = Ox)'(e"), a > 0, 
where gi are independent of e. Obviously, we obtain the relations 

which, by our assumption, imply 

gi^Q, z = l,2. 

Everything is different ^ if we assume that the coefficients g^, i = 1, 2, can 
depend on e. Here we consider only a special case of such dependence, which 
we shall use later. Namely, let 

5i = ^i + S'i(A(/./£), i = l,2, 

where Ai are independent of e and Si{p) decrease as \p\ — *■ oo. 
We assume that the estimate holds: 

\pSi[p)\ < const, i ~ 1,2. 

Let us find out what properties of the coefficients gi follow from the relation 

gi5{x - cj)i) + g2S{x - (/)2) = Ov'{e). 

Applying both sides of the equality to a test function ip, we obtain 

9if{(l)i) +g2f{(l)2) = 0(e) 

or, which is the same, 

[AM^i) + A^^ih)] + [SM(l>i) + S2V{<I>2)] = 0{e). (47) 

Let us consider the expression in the second brackets. Using Taylor's for- 
mula, we obtain 

[SMM+S2^i<l>2)] = SMM+S2^{(l>l)+S2{cl>2-M^'{^l+0(t>2), o<o<i. 

Now we see that 

52(A0/£)(02 - ^i) = {~pS2ip)}\^^^^/, ■ e = 0(£), 
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since the function pS2{p) is bounded uniformly in p G M'^. 
So we can rewrite relation (37) as 

AMM + A2ip{<j)2) + (Si + S2MM = Oie). 

Hence, as the coefficients Ai are independent of e, we, as usual, obtain 

Ai = 0, A2= 0, Si + 82^ 0. (48) 

Another method for analyzing relation (37) is the following. We assume that 
(j>i{t) are smooth functions, the relation (j)i{t*) — (p2{t*) holds for some t — t*, 
and, moreover, (f>'i{t*) ^ (t)'2{t*)- Then 

{Si5{x - (t>l),ip) + {S2S{X - (/)2), = Siip{x*) + S2(p{x*) 

+ SiO{t - 1*) + S20{t - 1*), X* = Mt*) = Mt*)- 

But 0{t - t*) - 0(A(/«). Therefore, we have 

(Ai + Si)6{x - cl>i) + {A2 + S2)6{x - 02) 

= Ai5{x - (t)i) + A25{x - </)i) + [Si + S2)5{x - x*) + Ov'{e). 

We again obtain relations (47). 



4.3 Complex germ lemma 



In this section, in the form convenient for us, we present the statement that 
plays an important role in Maslov's complex germ theory 00. 

Lemma 4.2. Let f{t) e , f{to) — 0, and f'{to) ^ 0. Let g{T, t) be a function 
that locally uniformly satisfies the estimates 



b'.9('''i^)l < const, \Tg'f.{T,t)\ < const, 
and g{T,to) — 0. Then the inequality 



-00 < T < 00, 



< Cts, 



where Ct — const, holds on any interval < t < T that does not contain zeros 
of the function f{t) except tQ. 

Proof: The fraction /(i)/(t — to) is locally bounded in t. The fraction r(7(T, i)/ 
(t — to) is also locally bounded. We have 
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t]{t-toy 



e ■ fit) ■ 



By the assumptions of the lemma, on the interval under study, the last multiplier 
on the right-hand side is bounded, while the product of the intermediate multi- 
plier and the expression in square brackets is bounded in view of the properties 
of the function ^(t, t). □ 
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Corollary 4.3. Suppose that the estimates in the assumptions of the lemma 
hold for < T < cx) (—00 < t < 0). Then the statement of the lemma 
holds on any interval [tQ,T] that does not contain zeros of the function f{t) and 
sgnr = sgn/(t), t e [to,T]. 

5 Justification of the weak asymptotic solution 

In this section we will prove that our weak asymptotic solution is in some sense 
"close" to the admissible weak solution of problem (1), (2). 

The existence of the admissible weak solution in our situation is obvious by 
Kruzhkov theorem (see Chapter 6). 

We will introduce admissibility conditions necessary for the uniqueness of 
the weak solution of considered problem. 

Definition 5.1. (Oleinik admissibility condition) We say that a weak solution 
u{t,x), t G M'^ , X e -K, of problem (1), (2) is admissible if it satisfies 

u+ = u{t, X* + 0) < u(t, X* — 0) — 

in every point of its discontinuity. 

Notice that such condition we can use only when the function u is piecevise 
continuous weak solutions of the considered problem for every fixed t € iR+. 
In that case Definition 15. II is equivalent to more general Kruzhkov admissibility 
condition (which can be applied on functions which are merely measurable): 

Definition 5.2. We say that the weak solution u{x,t), x E M, t E of 
problem (1), (2) is admissible if we have 

I I [dtipviu] + d^ipq{u)] dxdt + [ ij{x,0)r]{uoix))dx > 0, (49) 
Jo J m J R 

where q{u) — J r]'(u)f'{u)du and 77 e C'^{1R) is an arbitrary convex function. 

Using this definition, Kruzhkov proved the existancc uniqueness theorem 
(i.e. Theorem 6.2.2 in ^). 

We will prove that weak asymptotic solution tends in to the admissible 
weak solution of problem (1), (2). By the definition of the weak asymptotic 
solution for all G C°°([0, T]; C^(JRi)) we have: 

/ [u,t + ifius).]H^,t)dx = 0{e), 

JR 

uniformly in i G [0,r], T G 5?+. Integrating last expression with Jq dt we 
obtain: ^ 

/ / [uet + ifiu,)),](j)ix, t)dxdt = 0{e). (50) 

Jo J R 
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Now letting e — > we see that u(x,t) — w ~ lim Ue{x,t) is the weak solution 

of (1), (2). From the construction we see that u satisfy Oleinik admissibihty 
condition (since u is obviously piecevise continuous) and this implies Kruzhkov 
admissibility condition. Furthermore, it is easy to see that we have: 

/ [dtMue)+d^,Mue)]dxdt + [ ^{x,0)T]{uo{x))dx > eO{l), (51) 
Jn J Hi. 

where q(u) — J rj' {u)f' {u)du and rj G C^{1R) is an arbitrary convex function. 

Relation (50) holds by (49) and the smoothness of the function Ue{x,t) for 
e > 0. 

Now we can repeat the procedure from Theorem 6.2.2, page 87., to 
obtain: 

Theorem 5.3. Let and u satisfy (48) and (50), respectively. There exists 
s > depending only on [uq, U] (interval in which initial data take values) such 
that for every t G [0, T) and every r > we have: 

||w(-,t)-Ue(-,t)|Ui(i,i<,) < {r + st)-eO{l). 
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